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ABSTRACT 

We revisit the phenomenon of elements diffusion in the intergalactic medium (IGM) in 
clusters of galaxies. The diffusion is driven by gravity, concentration and temperature 
gradients. The latter cause thermal diffusion, which has been so far ignored in IGM 
studies. We consider the full problem based on the Burgers' equations and demonstrate 
that the temperature gradients present in clusters of galaxies may successfully compete 
with gravity, evacuating metals from cooler regions. Under the combined action of 
gravity and temperature gradients, complicated metallicity profiles with several peaks 
and depressions may be formed. For a typical cool core cluster, the thermal diffusion 
may significantly reduce and even reverse the gravitational sedimentation of metals, 
resulting in the depression in their abundance in the core. This may have implications 
for diagnostics of the low temperature plasma in the centers of clusters of galaxies. 

Key vifords: X-rays: galaxy clusters. 



1 INTRODUCTION 

Clusters of galaxies are the largest gravitationally-bound ob- 
jects in the Universe with thousands of galaxies residing in a 
dark matter-dominated potential. Hot intergalactic medium 
(IGM), constituting the majority of the baryonic matter in 
galaxy clusters, is a natural laboratory for plasma physics 
displaying a wide range of phenomena. Despite seemingly 
simple nature of galaxy clusters, understanding of physical 
processes in hot plasma of IGM is far from being complete. 

Of particular importance is the heavy elements distri- 
bution. Indeed, on the one hand it affects significantly the 
outcome of cosmological measurements via the mean molec- 
ular weight and thus the derived mass of the clusters (e.g. 
iQin fc WulbOOd: IChuzhov fc Nusseill2003l ') and the baryonic 
density fe.g. lEttori fc Fabianlbooi l On the other hand the 
X-ray spectrum of gas at low temperatures, T< 3 x 10^ K, is 
dominated by the emission lines, most important of which 
are associated with iron. The diagnostics of low tempera- 
ture plasma in clusters of galaxies is based primarily on line 
emission from heavy elements and the distribution of the 
latter may have direct implications for various p henomena, 
including the cooling flow problem (|Fabianlll994l ). 

Despite the importance of the heavy element distribu- 
tion, flat metallicity profiles are often assumed. This as- 
sumption may turn out to be too crude in some cases. In- 
deed, the starting point for the metallicity profile is the pri- 
mordial abundance. It is, however, quickly modified by a 
number of competing processes such as metals injection by 
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supernovae, mixing and diffusion. Although the effects of 
the latter process are rather controversial and its relevance 
to IGM in clusters of galaxies is still under debate, the sce- 
nario where the diffusion plays a significant role in forming 
the metallicity profile does not seem to be entirely unrealis- 
tic. 

The diffusion of heavy elements i n the intracluster 
medium has been considered by e.g. 'Fabian & PringQ 
(1977); Gilfan ov & Svunyaev (1984); Chuzhoy fc N usse^ 
(12003). Most studies focused on the gravitational separa- 
tion (gravitational sedimentation) and found that it may 
result in significant changes of the gas composition through- 
out the cluster. Given the presence of temperature gradi- 
ents in IGM, an additional type of diffusion may play a role, 
namely thermal diffusion. Although the importance of the 
latter in physics in general and in ste llar physics in particular 
has b een appreciated long ago (e.g. iThoul. Bahcall fc Loebl 
11994"), it was largely ignored in the case of ISM in galaxies 
and I GM in cluster of galaxie s. The exceptions are stud- 
ies bv IChuzhov fc Loebl Jiooi) who solved the full prob- 
lem, however c onsidered the case of isothermal cluster and 
IPeng fc Nagail (2008) who mentioned the potential relevance 
of the thermal diffusion effects in galaxy clusters. At first 
glance, simplified approach ignoring thermal diffusion may 
seem to be reasonable since the temperature variations in 
galaxy clusters are rather small. It turns out, however, that 
even moderate temperature gradient may be sufficient to 
compete the gravitational force. The consideration of the 
full problem including the effects of thermal diffusion is the 
subject of present paper. 

The paper is structured as follows. In section [2] we dis- 
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cuss the nature of the metalhcity profiles in galaxy clus- 
ters. In section [3] we describe the basic equations, numer- 
ical scheme and simulate the diffusion-driven evolution of 
a model galaxy cluster. Finally, in section [5] we discuss the 
obtained results. 



however, the time required to reach the equilibrium distri- 
bution may easily exceed the Hubble time and such a steep 
abindance profile is not formed. 

(iii) Concentration gradients. The diffusion driven by con- 
centration gradients works to restore the uniform abundance 
profile. 



2 METALLICITY PROFILES IN GALAXY 
CLUSTERS 

Although the observed metalhcity profiles are known to 
vary from cluster to cluster, they do exhibit som e sim- 
ilarities in their behaviour (e.g. IVikhlinin et al.l l2005l : 



[Sanderson. O'SuUivan fc PonmanlT2009l ). Typically, in cool 
core (CC) clusters the heavy element abundance rises from 
the cluster outskirts to the core. At the same time in the cen- 
tral parts, r < 10 — 20 kpc, a fraction of CC clusters show 
a depression in metalh city, whose nature a nd even reality is 
still under debate (e.g. IWerner et al]l2009l ). 

The metalhcity profiles in galaxy clusters form as a re- 
sult of interplay between different processes such as injection 
of metals from the supernovae-driven galactic winds, turbu- 
lent mixing and diffusion. The importance of the latter in 
plasma with turbulent magnetic field has been debated for 
a long time (e.g. iBakun in 2005; Lazarian 2007 ). As it is well 
known, a magnetic field hinders the diffusion across the field 
lines. However, given the chaotic character of the field in a 
turbulent medium, the efficiency of transport processes may 
be high enough to make the diffusion and t hermal conduc- 
tion important (jNaravan fc MedvedevllioOll ). In the follow- 
ing we assume that this is indeed the case. 

The diffusion is driven by the gravity, temperature and 
concentration gradients. Below we consider each in details. 

(i) Gravitational sedimentation. The equilibrium distri- 
bution of heavy elements, Us oc ^-"^sgh/kT ^ more concen- 
trated towards the cluster's center than that of the hydro- 
gen. Thus, due to the gravitational force, heavy elements 
tend to sink towards the cluster center, while light elements 
rise. 

(ii) Thermal diffusion. Thermal diffusion was discovered 
theor etically by Enskog and C hapman in the 1911-1917 
(e.g. IChapman fc Cowling! 1 19521 ') . A clear physical expla- 
nation of the effect has been given significantly later by 
iMonchick fc MasonI (|l967l '). who demonstrated that thermal 
diffusion results from the dependence of coUisional frequency 
on the velocity of particles. In the gas composed of two ele- 
ments, the effect of thermal diffusion is to evacuate the more 
massive and more highly charged species from the colder re- 
gion of the gas. In the case of a gas with three species - 
the abundant light ions, electrons and a small admixture of 
heavy ions with charge Zs and atomic masses As, the equi- 
librium distr ibution of the heavy element may be obtained 
analytically (|Burgerslll969l ): 



ns(r) cxr° 



exp 



(As - 0.5 Zs 



gjr) 

T{r) 



dr 



(1) 



where alpha is constant of the order of 3. This formula 
demonstrates the importance of thermal diffusion for heavy 
elements. Indeed, substituting Zs=26 for the iron, we obtain 
> 2000 for the power of T in the first multiplier. In practice. 



Another process affecting the heavy elements distri- 
bution is turbulent mixing which tends to wash out the 
metalhcity gradients. The effects of mixing depend on 
the spatial scales over which the turbulence is operat- 
ing and on the chaotic velocities spectrum. Mixing may 
also be driven by dynamical effects such as "sloshing" 
l|Ascasibar fc Marikevitchll2006l '). 

In the next section we construct a model cluster of 
galaxies and calculate the evolution of its metalhcity pro- 
file due to diffusion of elements. Our goal is to consider the 
full problem of diffusion for an idealized cluster in order 
to demonstrate the effects of heavy elements transport in 
galaxy clusters, with particular emphasis on the modifica- 
tions brought about by the thermal diffusion. Construction 
of a fully self-consistent model is beyond the scope of this 
paper, therefore we make several simplifying assumptions, 
namely: 

(i) We consider the idealized case of a static cluster as- 
suming that it's mass remains constant throughout the du- 
ration of simulations. In real case the cluster mass contin- 
iously grows via me rgers and accretion of DM halos (e.g. 
ICohn fc Whit3l2005l ). Since we are mostly interested in cool 
core clusters which are (most likely) not aff ected signifi- 
cantly by mergers in t he past few Gyrs (e.g. ISantos et al.l 
l2008l : |Poole" et al .1120081 ). this should not compromise our re- 
sults significantly. 

(ii) We assume that there is an exact balance be- 
tween heating and cooling of IGM which maintains the 
temperature profile constant throughout the simulations. 
This implies that the IGM heating via the thermal con- 
duction (which accompanies diffusion as these are the 
two sides of one phenomenon) and due to action of the 
central AGN equal the radi ative cooling. As shown by 
IGuo. Oh fc Ruszkowskil(|2008l ). a stable configuration of this 
kind may indeed exist. 

(iii) The IGM enrichment due to the continuous injection 
of metals from galaxies is ignored. 

(iv) The effect of mixing on the metalhcity profile of the 
cluster is ignored. 

The possible effect of the latter two processes is discussed in 
the section [431 



3 NUMERICAL EXPERIMENT 

We consider a spherically-symmetric cluster of galaxies with 
multicomponent plasma residing in a dark matter domi- 
nated gravitation al potential. The latter is assum ed to follow 
the NFW profile (|Navarro. Frank fc Whitil 19971 ). The clus- 
ter is assumed to be initially in the hydrostatic equilibrium. 
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3.1 Basic equations 

We consider AI species {M — 1 ions plus electrons) with 
concentrations ris, mean velocities Us, atomic masses As and 
charges Zs- All species have the same temperature Ts = T. 
The diffusion velocities Ws are defined as: 

Ws = Us — u, (2) 

where u is the mean fluid velocity: 

y^UsmsUs 
u= ^ . 

l^s 

The hydrodynamics of multicomponent gases is gov- 
erned by Burgers' equations (jBurgcrs 1969,). The Burgers' 
equations for the gas in the hydrostatic equilibrium (the 
mass, momentum and energy conservation equations) are 
(|Burgerslll969l : lThouL Bahcall fc Loeblll994l ): 



(3) 



dus 



+ 4r ^ [r'^ris (ws + w)] = 0, 



dr 



d{nskBT) 
dr 



+ UsTUsg — UsZseE - 



(4) 



Kstliwt ~ ws) + (5) 



+0.6{xstrs — Vstrt)], 



and 



5 , dT , 3 



-nsfcs— = > Kst{-Xst{ws - wt) - (6) 
2 dr ^ 2 

-yst\^&Xst{rs + rt) + Ystrs - 4.3a;strt]} - O.Sifssrs. 

Here g — [G'A/(r)/r^] is the gravitational acceleration mod- 
ulus, E is the radial electric field, Xst = ^st/rris, yst = 
fist/mt with ^J,st = msmt/{ms + nit) and Yst = 3yst + 
1.3xstmt/ms. Coefficients in equations correspond to pure 
Coulomb potential with a long-range cutoff at the De- 
bye length (|Thoul. Bahcall fc Loeblll994 ). The quantities 
are the s o-call ed residual heat flow vectors introduced by 
iBurgerd (|l969l ). playing the role of additional independent 
variables. The thermal diffusion arises due to the interfer- 
ence of equation [S] and the energy equation [S] via the rs . 
The friction coefficient between species s an t is: 

Kst = {2/S)fist{2kBT/fisty^^nsntcrst, (7) 

where the cross section between particles s and t is: 

ast = 2y^e*z!zUkBT)-^\nAst. (8) 

We assume the Coulomb logarithm to be InAat — 40 every- 
where. 

In addition, charge neutrality and the absence of cur- 
rents should be satisfied throughout the plasma. The former 
is not imposed explicitly in our code and is used to control 
the solution, while the latter is used to close the system of 
Burgers' equations (see below). 

3.2 Simulations 

The galaxy cluster is assumed to be initially in the hydro- 
static equilibrium, i.e. there is no net mass fiow, u = 0. 
However, once diffusion starts to operate, the total pres- 
sure changes and the gas is pushed out of the equilibrium, 
Vp + pg 7^ 0. Indeed, while heavy elements flow to the center 
of the cluster, hydrogen flows outwards and there is a net 



outflow of particles. The gas adjusts itself to the new equilib- 
rium configuration via the contraction on the hydrodynami- 
cal timescale. Given the much smaller timescale of the latter 
as compared to the diffusion timescale, the departure from 
equilibrium will be small enough for the hydrostatic equilib- 
rium form of Burgers' equations to be valid throughout the 
lifetime of a cluster (see also sec. I3.2.3|l . 

In order to solve the problem we use the following nu- 
merical scheme: 

(i) Based on the concentrations of elements nj"^ on k-1- 
th time step, solve Burgers' equations ^ and ((S} for the 
diffusion velocities Ws~^ ■ 

(ii) Given the diffusion velocities wj"^, the total mass 
ffow velocity u^~^ and ns~^, solve the continuity equation 
for the concentrations on the next time slice n*. 

(iii) Solve the Euler equation 



du 
It 



P 



(9) 



for the velocity u'' , restoring the hydrodynamic equilibrium 
in a cluster, using u*'~^ and n*"^. 

Below we discuss each step in detail. 



3.2.1 Diffusion velocity calculation. 

Burgers' equations (eqs [5] and IB} present a linear inverse 
problem of the form Axy = b on the diffusion velocities Wa , 
heat vectors of species r^ and electric fleld E. The diffusion is 
driven by the gradients VT, Vris, Vp(g(r)), with all these 
quantities residing on the r.h.s. b. Given the M species we 
have 2M-I-1 unknowns and thus need 2M-I-1 independent 
equations. Burgers' equations give us only 2M-1 independent 
equations, because diffusion velocities appear in eqs. [5] and 
[B]only in the form {ws — Wt). The two additional conditions 
needed to close the system of Burger's equations are defined 
as follows: 



E 



piWi 



0, 



(10) 



which follows from the definition of the diffusion velocity 
(eqs[2]and Finally the absence of electric currents gives: 



ZiUiUi = 



(11) 



The resulting system of equations is solved for the 
Ws and rs using the sing ular value decomposition (SVD) 
method (|Press et al.lll992 



3.2.2 Continuity and Euler equations. 

The problem we are solving is a mixed advection-diffusion 
problem. We have found that it can be solved with sufficient 
accuracy using the combination of forward- time centered- 
space (FTCS) and Lax methods (|Torolll999l '). with the for- 
mer applied to the continuity equation and the latter to 
Euler equation. The Lax method is known to introduce nu- 
merical diffusion; however, we performed a number of tests 
and found that it is not crucial for our case (see Sec. I3.2.3|) . 

We use a homogenous grid spanning from to 2 Mpc 
with step dr = 0.25 kpc. The timestep is chosen to be small 
enough to ensure the stability of the scheme, dt — 2 — 5 ■ 
10"^ Gyr. 
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Since we are mainly interested in the diffusion in the 
central region of the cluster, we are relatively free in the 
choice of the outer boundary condition. We assume that the 
galaxy cluster is immersed in the medium with constant den- 
sity, i.e. there is an infinite reservoir of gas which maintains 
the gas density together with the element mass fractions at 
the outer boundary constant. In calculating velocities at the 
outer boundary we use one-sided derivarives. 

At the inner boundary r = the continuity equation 
has a singular point. Therefore we expand eq. Q in the 
vicinity of zero and obtain the folowing expression for the 
concentration at the k-th time step: 



k k-1 



3 X {dt/dr)n'l ..'^{u2 ^ +-u^2,/) 



(12) 



We also required that the density gradient should be equal to 
zero at r = and used a one-sided derivative for the velocity 
gradient Vit. In addition all velocities should be zero at the 
inner point: Ws = and u = 0. 



3.2.3 Code validation. 

We further performed a number of tests to check the validity 
of our code. 

(i) Numerical diffusion. We obtained the solution replac- 
ing the Lax method with the FTCS. The resulting code 
is free from numerical diffusion, however the FTCS-based 
scheme is less suited for our purpose due to a well-known 
instability of FTCS for hydrodynamical problems. On the 
other hand since the hydrodynamical terms in our equations 
appear only as a result of a slow diffusion process, we were 
able to obtain a stable solution for the 7 Gyr simulation. The 
solutions obtained using the combination of Lax-|-FTCS and 
FTCS-I-FTCS schemes are consistent with each other within 
a few percent. 

(ii) Boundary conditions. We check that the solution is 
insensitive to the boundary conditions. For this purpose we 
run our code for 7 Gyr with different outer radii and found 
no significant difference. We also considered the opaque 
outer boundary, setting to zero all the velocities at the outer 
radius. The solution at t = 7 Gyr in this case is also con- 
sistent with the one for the default boundary conditions al- 
most everywhere, within a fraction of a percent. Solutions 
obtained with different boundary conditions differ from each 
other only within ~ 100 kpc from the outer boundary. This 
dependence is unimportant and can be neglected. 

(iii) Hydrostatic equilibrium. We also checked that the 
departures from the hydrostatic equilibrium are not too 
large, for the Burgers' equations to be valid, namely that 
IVp/p + gl <C IVps/ps +g\ throughout the duration of simu- 
lations. This condition is satisfied with an accuracy of better 
than a fraction of a percent. 

(iv) Electro-neutrality. Since we do not use the elec- 
troneutrality condition explicitly in our code, we check 
that the final solution satisfies it. We obtain that 
\J2Z^n,/n, - l| < IQ-^ 

(v) Comparison with an analytical solution. In order to 
check the implementation of the Burgers' equations in the 
code, we consider a plasma composed of protons, electrons 
and a small admixture of iron in the hydrostatical equilib- 
rium. Solving the Burgers' equations analytically we obtain 
the coefficient a in the power of temperature in eq. [1] a ~ 



3.1, and compute the equilibrium distribution of iron. For 
this distribution we calculate the diffusion velocities using 
our code and verify that they equal zero as expected for an 
equilibrium distribution. We obtain, that wfb ^ 10"'^ km/s. 



4 RESULTS AND DISCUSSION 
4.1 Model cluster 

For the temperature profile of the cluste r, T(r), we use 
the analytical best-fit model for A262 from IVikhlinin et al.l 
( 2006. ). Outside the range of its applicability the profile was 
extended as illustrated in the Figure [1] The total mass den- 
sity of the cluster is assumed to follow the NFW profile 
iNavarro. Frank fc Whit3l 19971 '): 



pNFW ~ 



(r/r,)(l+r/r,)2- 



(13) 



We tak e rs — 184 kpc, the observed value for the A262 
cluster IVikhUnin et aLlbOOd ). 

Given the T{r) and rs and assuming hydrostatic equi- 
librium, we determine the value of ps resulting in gas density 
pg{r) which follows as close as possible the analyt ical best-fit 
model for gas density from lVikhlinin et al.l (|20061 ). With the 
obtained value of normalisation, ps = 1.03 x 10~^® g/cm^, 
the total mass inside the radius r^ = 70 kpc obtained inte- 
grating the pNFw(r) equals the total mass calculated using 
the temperatu re and gas densi t y pro files (and their extrap- 
olations) from IVikhlinin et al.l (120061 ). assuming the hydro- 
static equilibrium. The overall agreement between these two 
masses is better than ~ 20 — 30% in the range of 20-450 kpc 
with the limits reflectin g the bounds of the re gion with a suf- 
ficient quality of data (|Vikhlinin et al.ll2006l ). The final gas 
density profile is calculated using the adjusted NFW profile 
pNFwir), temperature profile T(r) and assuming the hy- 
drostatic equilibrium. The normalization of the gas density 
profile corresponds to a gas mass fraction within rsoo (the 
radius where the overdensity is equal to 500) of /500 = 0.14. 
The resulting total mass, gas density and temperature pro- 
files are shown in fig. [T] 

To illustrate the role of the thermal diffusion, we also 
consider an isothermal cluster with temperature kT = 2 
keV and the same total mass density profile. The initial gas 
density profile is calculated similarly to above, assuming the 
hydrostatic equilibrium (fig. 

We consider a 4-component plasma consisting of 
He"'"^, Fe"""^^ and electrons. The initial metallicity profiles are 
assumed to be flat with the mass fractions of hydrogen and 
iron being 0.75 and 1.8 • 10"'^ respectively. To demonstrate 
the effects of diffusion for different species, we also perform 
calculations replacing the Fe^^'' by O^*. 



4.2 Diffusion in clusters of galaxies: temperature 
vs. gravity 

We plot on fig. [1] and [2] the diffusion velocity profiles at the 
initial moment of time for a cool core and isothermal models 
respectively. Positive values correspond to outfiow of parti- 
cles, negative to inflow. As is clear from the plots, the pres- 
ence of temperature gradients drastically changes the char- 
acter of diffusion. While in an isothermal cluster the effect 
of diffusion is restricted to sedimentation of He and heavy 
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Figure 1. Left: Dark Matter and gas density profiles and temperature profile of the model galaxy cluster. Right: Corresponding initial 
diffusion velocity profiles for H (red), He (green), Fe (blue) and O (magenta, dashed line). Positive velocities correspond to outflow of 
particles, negative — to inflow. A positive peak at small radii and a negative one at the cluster outskirts are caused by the thermal 
diffusion and correspond to strong temperature gradients in the IGM (see the left panel). 




nuclei, in the non-isothermal case the picture is more com- 
plicated. The combined effect of gravity and temperature 
gradients leads to formation of a complex diffusion veloc- 
ity profile with prominent features at the positions of strong 
temperature gradients (~20 kpc and ~ 400 kpc in fig.[T|. At 
small radii, the sharp decrease of temperature counteracts 
gravity and reverses the diffusion flow, evacuating heavy el- 
ements from the cool core of the cluster. On the other hand, 
the negative temperature gradient in the outer part of the 
cluster works in the same direction as the force of gravity, 
increasing the velocity of gravitational sedimentation. Note 



that due to the lower temperatures and higher gas density, 
the diffusion velocities in the centers of cool core clusters 
are smaller than in higher temperature IGM usually consid- 
ered in the context of the gravitat ional sedimentation of ele - 
ments in clusters of galaxies (e.g. IChuzhov fc NusserlboOSi ). 
The goal of the present analysis is to study the possible role 
of thermal diffusion in the IGM. We therefore focused on 
the cool core clusters having strong temperature gradients, 
rather than on the high temperature clusters characterized 
by higher diffusion velocities. Because of the relatively small 
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Figure 4. Same as in Fig. E] but for an isothermal cluster from Fig. [2] 



diffusion velocity, ions can drift only through a small fraction 
of the cluster, < 20 kpc, during it's lifetime. 

The overall effect of diffusion is further illustrated by 
the metallicity profiles after 3 and 7 Gyrs shown in figure |3l 
At large radii the heavy elements flow inwards due to com- 
bined effect of thermal diffusion and gravity. The tempera- 
ture peak at ~ 300 kpc leads to the enhancement of heavy 
elements abundance at this radius. At intermediate radii the 
temperature gradient is relatively weak and thus metals flow 
under the action of gravity to the center. Due to the strong 
temperature gradient at r ~ 20 kpc, the inward diffusion 
of heavy elements stops and the second peak is formed in 
heavy elements distribution. At the same time metals from 



the core of the cluster ffow up the temperature gradient, 
making the peak more prominent and creating a depression 
in the metallicity at smaller radii. Because of the nearly flat 
temperature distribution in the inner ~ 10 kpc, assumed 
here, the abundance at the center of the cluster is also some- 
what enhanced by the gravitational sedimentation, with the 
effect being stronger for helium, than for iron. In contrast 
to this, the picture is much simpler in the isothermal cluster 
model shown in flg. [5] In this case the density proflle evo- 
lution is driven by the gravitational sedimentation of heavy 
elements leading to a monotonic increase of their abundance 
towards the cluster center. 

In figure [5] we compare the diffusion velocity profiles 
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Figure 5. Diffusion velocity profiles for Fe (blue line) and He 
(green) in the model galaxy cluster after 7 Gyr. Dashed lines show 
the initial diffusion velocities. 

after 7 Gyr with the initial distributions. Clearly, the outflow 
of Fe at small radii is only shghtly reduced, while for He it 
almost disappears. The diffusion velocities outside the inner 
parts of the cluster remain almost unchanged. 

4.3 Dependence on the atomic mass and charge 

Both gravitational sedimentation and thermal diffusion ve- 
locities depend on the ion species. As is clear from Fig. (2] 
the gravitational sedimentation in an isothermal cluster is 
the most effective for the helium ions, with the higher Z 
ions such as Fe and O diffusing somewhat slower. At the 
same time, the diffusi on velocities are almost e qual for the 
latter two, wpe ~ wo l|Chuzhov fc Nusseill2003l ). For a non- 
isothermal temperature profile, the equilibrium distribution 
of the heavy element depends on the element's charge with 
more highly charged species being more concentrated to- 
wards peaks in temperature profile of the cluster (e.g. eq.[l|. 
As a result, the abundance profile of He appears to be 
smoother than that of Fe (fig. On the other hand, on 
the timescales of interest, the effect of diffusion for elements 
heavier than He is identical from the observational point 
of view. The strong charge dependence of thermal diffu- 
sion (eq. [ij becomes relevant only on significantly longer 
timescales, due to the small VT in the 30-200 kpc interval. 

4.4 Dependence on the cluster parameters 

The overaU pattern of diffusion depends critically on the 
cluster parameters, due to the strong temperature, density 
and gravity dependence of diffusion velocities. Since the typ- 
ical temperatures in clusters of galaxies vary from less than 
a 1 keV up to ~ 10 keV and the diffusion velocity scales 
with temperature as w on T^^' (e.g. lThoul. Bahcall fc Lo^ 



1 1994 ). the variations in the diffusion velocity may exceed an 
order of magnitude. On the other hand, the temperature, gas 
density and gravi ty in clusters of galaxies are by no means 
independent (e.g. IVikhlinin et al.|[200Q ) . Thus the character 
of the diffusion is a result of the complex interplay between 
the cluster parameters. 

The reversal of the diffusion fiow for He and metals 
at small radii (fig. [1} is a consequence of a large gradient 
in the temperature profile, which counteracts the force of 
gravity. The balance between these two forces is sensitive 
to the details of the temperature distribution. Rather small 
and otherwise almost insignificant modifications of the tem- 
perature profile may lead to the diffusion pattern different 
from that shown in fig. [1] and [S] To illustrate this, we con- 
sider a cluster of galaxies with an (ad hoc) smoother tem- 
perature profile, characterized by smaUer gradients in the 
cluster core, but with the same total mass distribution. The 
cluster parameters are shown in fig. |6] As is obvious from 
the figure, the initial distribution of the gas density does 
not differ significantly from our default model. However the 
diffusion velocity profile changes dramatically - the positive 
velocity peaks disappear for both elements and He velocity 
becomes negative in the entire range of radii (fig. [6]). The 
resulting metallicity profile is correspondingly significantly 
smoother (cf. fig.|3]). However, thermal diffusion counteracts 
the effect of the gravitational sedimentation in this case as 
well, significantly reducing its effect for He and resulting in 
a remarkable central depression in the Fe abundance. 

Further decrease of the temperature gradient in the 
cluster core can make gravity the dominant factor. Based 
on the solution of Burgers' equations we estimate that for 
the gravitational potential considered here, the thermal dif- 
fusion equals to the force of gravity for the temperature 
gradient of AT ~ 0.1 keV per 20 kpc. Note that this is ~2 
times smaUer than the temperature gradient in the profile 
shown in the fig. AT ~ 0.2 keV per 20 kpc at r = 20 kpc. 

4.5 Turbulent mixing and metals injection from 
cD galaxy 

An idealized model constructed in this paper ignores turbu- 
lent mixing and metal enrichment of the IGM by the central 
cD galaxy - the two important factors which potentially may 
have a significant impact on the metallicity profiles. The 
turbulent mixing tends to smoothen the abundance distri- 
butions, while the cD galaxy increases the heavy element 
abundances in the inner parts of the cluster due to metals 
produced by (mostly) SNIa explosions. These two factors 
are discussed below. 

In order to estimate role of the cD galaxy, we 
measure the K-band luminosity of A262 using the 
image from 2MASS Extended Source Image Server 
(http:/ /irsa. ipac.caltech.edu / applications / 2 M ASS / Pub- 
GalPS/) and based on the calibration of iMannucci et al] 
(,2005< ) for E/SO galaxies convert it to the SN la rate. 
Combining the latter with the iron yield of rj ~ 0.7 Mq 
l|lwamoto et al.|[l999l ) and assuming perfect mixing of the 
SNIa ejecta with IGM, we obtain the iron injection rate 
of 8.0 • 10"^ Me/yr within the inner ^ 19 kpc of the 
cluster. On the other hand, given the diffusion velocity, 
wfc ~ 10"" cm/s, and concentration of iron at this radius 
(see fig. [1]), we obtain the diffusion-driven iron outfiow 
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Figure 6. Left: Dark Matter, gas density and temperature profiles in the model galaxy clusters witli an (ad hoc) temperature distribution 
with a moderate gradient in the core. The total mass distribution is the same as in our default model shown in fig. [l] Right: Diffusion 
velocity profiles for He (red) and Fe (blue) at the initial moment of time. For comparison, the dashed lines in both panels show profiles 
for the case of a large temperature gradient from figH] Unlike the latter case. He fiows inwards over all radii. 
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Figure 7. Left: Abundance profiles for Fe (blue solid line), and He (green line) in the model galaxy cluster from fig.|6]after 3 Gyr. Right: 
The same, but after 7 Gyr. 



from the same region of ~ 2.3 ■ 10~^ M0/yr, which is 
several times smaller than the iron injection rate. The 
relation between the outflow and injection rates is opposite 
for a-elements. Indeed, assuming the oxygen yield of 
~ 0.1 M0 and solar metallicity, the oxygen injection rate is 
~ 1.1 ■ 10~^ M0/yr, which is several times smaller than the 
outflow rate, ~ 4.9 • 10~^ Me/yr. Taking these numbers at 
the face value, we may conclude that the impact of the cD 
galaxy on metallicity profiles varies from significant for iron 
peak elements to negligible for a-elements. 



We assumed above that complete mixing of the Super- 
nova ejecta with IGM takes place. However, the details of 
the SN ejecta interaction with the ICM are poorly under- 
stood. Indeed, the assumption of the total mixing leads to 
unreasonably high iron abundance in the cluster core, as il- 
lustrated by the following calculation. Assuming that the 
iron injection rate in A262 follows the stellar mass profile 
and taking into account the metal injection by Supernovae 
via an additional source term in eq.|4l we repeated the calcu- 



lation of the metallicity profile evolution described in sec. [3rl 
This resulted in a strongly peaked iron profile (Fig[8]) with 
unrealistically high metaUicity in the center of the cluster, 
Z ~ 7 Zq within the inner 10 kpc. Such a high iron abun- 
dance in the cluster core contradicts to observations of clus- 
ters of galaxi es, in particular t o the iron distribution in A262 
measured bv lVikhlinin et al] (|2005h . pointing at the need of 
a more sophistica ted model. Similar results were obtained 
bv lRebusco et ahl ((2006), who found the observed metallic- 
ity profiles in galaxy clusters to be broader than predicted 
under the assumption that the rate of metals injection fol- 
lows stellar light distribution. A possible explanation of such 
a discrepancy is offered by the turbulent mixing driven by 
the central AGN activity. This requires turbulent diffusion 
coefficients of 6 ■ 10^'' - 5 ■ 10^^ cm /s, corresponding to mixing 
length scales of 1 — 35 kpc and gas velocities « 200—400 km/s 
iRebusco et al.l [20061 ). However, the role of this process in 
the particular case of A262 and other cool core clusters 
is unclear. Indeed, th e presence of steep temper ature and 
abundance gradients l|Vikhlinin et al.ll2005l , 120061 ) rules out 
any significant mixing operating on spatial scales exceed- 
ing ~ 10 kpc, thus excluding a large part of the parameters 
space required for the turbulent mixing interpretation. 

To conclude, the role of mixing and metal enrichments 
remains to be controversial. An attempt to include them in 
the model in a straightforward manner leads to an appar- 
ent conflict with observations. This suggests that the role of 
these processes may be overestimated. Secondly, this points 
at the need of more sophisticated models. Consideration of 
such models will be a subject of the forthcoming paper. 



5 CONCLUSIONS 

We considered the effects of element diffusion in the inter- 
galactic medium on the metallicity profiles in clusters of 
galaxies. The diffusion processes in the IGM are driven by 
the gravity, concentration and temperature gradients. Grav- 
ity forces metals to sink towards the center, concentration 
gradients act to restore the uniform distribution, and the 
temperature gradients tend to evacuate heavy elements from 
the low temperature regions in the IGM. So far, attention 
has been paid to gravitational sedimentation, while thermal 
diffusion has been largely excluded from consideration. We 
demonstrate that the latter may drastically change the pat- 
tern of the diffusion and thus may play a significant role in 
formation of the metallicity distribution in IGM. 

We performed numerical simulations of the diffusion in 
the IGM of clusters of galaxies based on the full set of Burg- 
ers' equations. In order to illustrate the effect of thermal 
diffusion we considered three different models of cluster of 
galaxies with the same total mass profile but different tem- 
perature distributions: isothermal, and with large and mod- 
erate temperature gradient (figs. [1] (2] |6]). The evolution of 
metallicity profiles are shown on figs. [S] |4] and [7] As illus- 
trated by these plots, the combined effects of the tempera- 
ture gradients and gravity may result in a rather complex 
metallicity profiles in clusters of galaxies. 

1 Note, that due to the limited coverage of A262 by 2MASS, the 
metal injection rate at r > 19 kpc is somewhat underestimated. 
This, however, does not change our conclusions. 
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Figure 8. Abundance profile for Fe in the model galaxy cluster 
with injection of metals by cD galaxy after 7 Gyr. See section[43] 
for details. 



The efficiency of transport processes in IGM is still un- 
der debate. However, if element diffusion does operate in 
clusters of galaxies, its overall effect on metallicity distribu- 
tions may be opposite to what was previously thought, based 
on simplified considerations of isothermal IGM. Our results 
show that even moderate temperature gradients, ~ 0.5 keV 
per 100 kpc may be sufficient to counteract the force of 
gravity and to stop and even to reverse the diffusion flow 
of metals, thus preventing enhancement of heavy element 
abundances in the center due to gravitational sedimenta- 
tion. It is especially relevant to the cool core clusters, where 
the thermal diffusion may lead to depletion of metals in the 
cluster core. The implications of this process for the diag- 
nostics of the low temperature plasma in the centers of clus- 
ters of galaxies will be considered in the forthcoming paper 
(Shtykovskiy & Gilfanov, in preparation). 
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